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Introduction 

Knowledge  of  wave  and  surf  conditions  can  play  a 
major  role  in  planning  and  executing  a  successful 
amphibious  assault  on  a  targeted  beach  (Figure  1).  The 
dynamic  surf  zone  is  affected  by  waves  that  may  have 
propagated  from  long  distances  (swell)  and  by  refrac¬ 
tion,  which  causes  waves  to  become  aligned  with 
depth  contours  as  they  approach  the  coast.  Local  wind 
waves  and  tides  can  also  change  the  character  and 
width  of  the  surf  zone.  Local  bathymetry  (e.g.  beach 
profiles)  can  change  dramatically  due  to  wave  action 
from  storms.  Offshore  bars  can  form  and  migrate  due 
to  seasonal  changes  in  wave  dynamics  and  surf  condi¬ 
tions.  Accurate  predictions  of  surf  processes  are 
dependent  on  models  that  contain  the  appropriate 
physics,  and  realistic  environmental  inputs  to  these 
models. 

A  wave-tide-surf  modeling  system  has  been  devel¬ 
oped  to  support  naval  operations  in  littoral  waters.  The 
Integrated  Ocean  Prediction  System  (IOPS)  provides 
integrated  wave  information  from  deep  to  shallow 
water  into  the  surf  zone.  IOPS  is  composed  of  four 
component  models:  (1)  the  Wave  Action  Model 
(WAM);  (2)  two  shallow-water  wave  models  (steady- 
state  wave  model,  STWAVE,  and  the 
REFraction/DIFfraction  Model,  REFDIF);  (3)  two  tidal 
models  (ADvanced  CIRCulation  Model,  ADCIRC,  and 
the  US  Navy  Relocatable  Tide /Surge  Prediction 
System-PCTIDES);  and  (4)  the  Navy  Standard  Surf 
Model  (NSSM).  Different  physical  processes  and  spa¬ 
tial  scales  affecting  wave  conditions  across  shallow- 
water  regions  are  considered  when  these  models  are 
used.  In  addition,  the  modeling  techniques  delineate 
surf  characteristics  that  may  vary  substantially  along 
stretches  of  coast. 

System  Overview 

IOPS  is  a  component  of  the  Rapid  Ocean  Analysis 
Modeling  Evaluation  Relocatable  (ROAMER)  System, 
that  is  designed  to  provide  the  Naval  Oceanographic 
Office's  (NAVOCEANO's)  Warfighter  Support  Center 
(WSC),  with  the  capability  to  set  up,  run,  and  monitor 
model  performance  using  a  Graphical  User  Interface 
(GUI).  ROAMER  is  composed  of  a  relational  database 


in  which  information  (grids,  etc.)  about  each  model  is 
stored.  The  ROAMER  system  provides  the  user  access 
to  NAVOCEANO's  DBDBV  bathymetry  (global  5 
minute  resolution)  and  a  host  of  Fleet  Numerical 
Meteorology  and  Oceanography  Center  METOC  prod¬ 
ucts.  ROAMER  is  used  to  pull  in  the  required  input 
data  and  boundary  conditions  and  build  a  script  to 
launch  the  appropriate  model  run.  The  NAVOCEANO 
WSC  utilizes  the  NAVOCEANO  Major  Shared 
Resource  Center  (MSRC),  a  network  of  computers  that 
can  support  parallel  codes  and  multi-processors. 

The  IOPS  model  suite  is  depicted  in  Figure  2.  The 
NSSM  requires  four  inputs:  1)  shallow- water  wave 
spectra,  2)  water  elevation  from  a  tidal  model,  3)  beach 
profile(s)  and  4)  10-m  surface  wind.  The  regional  WAM 
provides  directional  wave  spectra  that  are  used  as  an 
offshore  boundary  condition  for  a  shallow-water  wave 


Figure  1.  Assault  Amphibian  Vehicles  from  AAV 
Platoon ,  Battalion  Landing  Team  3/1,  13th  Marine 
Expeditionary  Unit  (Special  Operations  Capable),  made 
the  most  of  their  visit  to  Iwo  Jima  by  landing  all  AAVs 
onto  Green  Beach,  Jan.  26,  2001.  Marine  photo  by  Sgt. 
M.C.  Brown. 
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Figure  2.  IOPS  model  components.  The  Navy  Standard  Surf  Model  requires  four  inputs:  directional  wave  spectra  from  a  shal¬ 
low-water  wave  model ,  water  elevation  from  a  tidal  prediction  model ,  beach  profiles ,  and  wind  speed  and  direction. 


model.  A  tidal  prediction  model,  which  can  be 
forced  with  high-resolution  winds,  provides  tidal 
inputs  to  NSSM.  Last,  beach  profiles  and  winds  are 
passed  to  NSSM. 

WAM 

WAM  (WAMDI  group  1988;  Komen  et  al.,  1994) 
was  originally  developed  at  the  Max  Planck  Institute 
for  Meteorology  in  Hamburg,  Germany,  and  updated 
by  a  consortium  of  collaborating  scientists  in  Europe  to 
the  current  version  known  as  WAM  Cycle  4.  The  U.S. 
Army  Corps  of  Engineers  Waterways  Experiment 
Station  (WES)  at  the  Coastal  Hydraulics  Lab  (CHL)  in 
Vicksburg,  Mississippi,  adapted  command  files  and 
source  code  for  use  by  the  U.S.  Navy.  At  Vicksburg, 
efforts  continue  to  modify  WAM  (Jensen  et  al.,  this 
issue)  for  more  efficient  parallel  computing. 

NAVOCEANO  made  further  modifications  to 
adapt  the  model  for  efficient  operational  processing  on 
the  NAVOCEANO  MSRC.  The  primary  computational 


unit  for  the  numerical  model  is  the  Cray  Scalable- 
Vector  (SV1)  platform,  which  can  perform  vector  calcu¬ 
lations  while  changing  the  number  of  processors  on  the 
fly.  While  WAM  is  coded  for  vector  processing,  the  scal¬ 
able  aspect  of  the  SV1  allows  for  running  the  model  on 
as  many  as  eight  CPUs.  As  a  result,  NAVOCEANO  can 
run  a  48-hour  WAM  forecast  for  the  entire  globe  in  less 
than  30  minutes  and  has  been  able  to  run  45  separate 
domains  around  the  world  every  12  hours. 

In  addition  to  high  throughput,  it  is  possible  to  run 
WAM  in  a  large  variety  of  domains  of  differing  shapes 
and  sizes,  due  to  modifications  that  allow  for  dynamic 
memory  allocation.  Using  Fortran90,  a  single  exe¬ 
cutable  program  can  be  run  for  any  size  without  chang¬ 
ing  the  source  code  and  recompiling.  In  addition,  rather 
than  having  to  run  a  separate,  individual  host  model  to 
provide  boundary  conditions  for  each  nest  (Figure  3a), 
the  code  allows  for  a  single  host  to  output  boundary 
conditions  for  multiple  nests  (Figure  3b),  which  can  be 
changed  on  the  fly.  In  this  way  model  domain  setup  can 
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Figure  3.  WAM  nesting  a)  stovepipe  approach  to  model  domain  set  up  in  which  a  host  is  required  for  every  nest,  b)  improve¬ 
ments  allow  multiple  nests  to  he  run  from  a  single  host. 


be  accomplished  in  a  few  minutes.  Also,  more  and 
larger  domains  can  run  on  the  computational  machine 
in  less  time. 

For  Navy  applications.  Basin  scale  WAM  runs  are 
forced  with  the  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS),  which  has  a  horizontal 
resolution  of  1.0°.  For  embedded  nests  of  higher  reso¬ 
lution  (generally  <  A0)  winds  from  the  Coupled  Ocean 
Atmospheric  Mesoscale  Prediction  System  (COAMPS) 
are  used  when  available. 

Regional  domains  cover  large  areas  such  as  the 
entire  Mediterranean  Sea  at  a  0.1°  grid  spacing  and  as 
small  as  1-minute  domains  about  the  Hawaiian 
Islands.  Other  regions  include  the  seas  around  North 
and  Central  America,  the  Indian  Ocean,  the  Asian 
Pacific,  Indonesia,  and  seas  surrounding  the  Australian 
continent.  WAM  output  is  then  post-processed  to  pro¬ 
vide  over  a  hundred  sets  of  tailored  graphics  products 
on  the  World  Wide  Web  (Figure  4,  see  also 
http://www.navo.navy.mil)  for  access  to  the  warfight¬ 
er.  Products  include  a  rendering  of  significant  wave 
height,  primary  wave  direction,  swell  height  and  direc¬ 
tion,  and  wave  period.  These  parameters  are  calculated 
from  the  model's  computational  arrays'  (of  discrete  fre¬ 
quency  and  direction)  spectral  values.  These  values  are 
passed  from  host  to  nests  as  boundary  conditions. 
Also,  WAM  passes  spectral  values  as  boundary  condi¬ 
tions  to  shallow-water  wave  models  such  as  ST  WAVE. 

WAM  at  NAVOCEANO  has  been  enhanced  for 
increased  throughput  and  greater  flexibility  to  run  on 
many  domains  and  to  provide  the  input  for  littoral, 
shallow-water  model  predictions  anywhere  in  the 
world.  In  the  future,  WAM  will  be  ported  on  to 
new  architectures  such  as  the  IBM-SP  and  LINUX  PCs. 
These  increased  capabilities,  current  and  planned. 


will  be  integrated  into  the  IOPS  wave-tide-surf 
modeling  system. 

Shallow-Water  Wave  Models 

The  IOPS  model  suite  utilizes  two  shallow-water 
wave  model  options:  (1)  STWAVE  and  (2)  REFDIF. 
STWAVE  can  generate  2-D  outputs  of  significant  wave 
height,  peak  period,  and  mean  wave  direction. 
Directional  wave  spectra  saved  at  user-specified  loca¬ 
tions  can  be  used  for  generating  subsequent  surf  fore¬ 
casts.  In  another  IOPS  application,  REFDIF  can  gener¬ 
ate  shallow-water  wave  spectra  when  wind  is  not  a 
factor,  since  the  model  does  not  include  wind  effects. 
Utilizing  a  linear  transfer  coefficient  technique,  WAM 
deep-water  wave  spectra  are  transformed  to  shallow- 
water  wave  spectra.  While  this  approach  can  require 
extensive  computations  to  derive  these  coefficients,  the 
payoff  is  at  operational  run-time,  when  accurate  results 
can  be  produced  in  a  few  minutes,  versus  hours.  The 
following  section  gives  a  more  detailed  description  of 
these  two  shallow-water  wave  models. 

STWAVE 

The  Steady-state  spectral  Wave  model  (Resio  1987, 
1988a,  1988b;  Smith  et  al.,  1999)  is  a  steady-state  finite- 
difference  model  based  on  the  wave-action  balance 
equation.  STWAVE  is  a  half-plane  model,  which  con¬ 
siders  wave  energy  propagating  from  the  offshore 
toward  the  nearshore  from  the  x-axis  of  the  grid,  which 
is  typically  the  approximate  shore-normal  direction. 
STWAVE  can  simulate  depth-induced  wave  refraction 
and  shoaling,  current-induced  refraction  and  shoaling, 
depth-  and  steepness-induced  wave  breaking,  diffrac¬ 
tion,  wind-induced  wave  growth,  and  wave-wave 
interaction  and  white  capping  that  redistribute  and  dis- 
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Figure  4.  Standard  ocean  surface  wave  forecast  products  used  by  the  warfighter  as  provided  by  the  Warfighting  Support 
Center  at  the  U.S.  Naval  Oceanographic  Office. 


sipate  energy  in  a  growing  wave  field. 

STWAVE  assumes  that  offshore  wave  conditions 
are  homogenous.  Thus,  the  wave  energy  (provided 
from  a  regional  WAM)  is  held  constant  on  the  domain's 
offshore  boundary.  The  steady-state  formulation 
reduces  computational  time  and  is  appropriate  for 
wave  conditions  that  vary  more  slowly  than  the  time  it 
takes  for  waves  to  travel  across  the  computational  grid. 

STWAVE  includes  a  wave-current  interaction  fea¬ 
ture  based  on  a  barotropic  or  depth-averaged  current. 
If  strong  vertical  gradients  in  the  current  occur,  their 
modification  of  shoaling  and  refraction  will  not  be  rep¬ 
resented  in  the  model. 

A  finite-difference  model,  STWAVE,  is  formulated 
for  a  Cartesian  coordinate  system  using  square  grid 
cells  (Ax  =  Ay).  STWAVE  operates  in  a  local  coordinate 
system  with  the  x-axis  oriented  in  the  cross-shore 
direction  and  the  y-axis  oriented  along-shore.  The  y- 
axis  is  usually  aligned  with  the  bottom  contours  along 
the  outer  boundary  of  the  grid.  In  a  standard  IOPS 


implementation,  a  STWAVE  grid  may  be  set  up  for  an 
area  of  coastline  ranging  from  15-30  km  along-shore 
and  10-20  km  in  the  cross-shore.  To  run  a  subsequent 
surf  model  application,  a  series  of  locations  is  specified 
in  which  shallow- water  wave  spectra  are  to  be  saved. 

REFDIF 

The  REFDIF  model  of  Kirby  and  Dalrymple  (1994) 
incorporates  the  effects  of  shoaling,  refraction,  energy 
dissipation,  and  diffraction  on  the  propagation  of  water 
waves  over  a  field  of  irregular  depths  and  around  irreg¬ 
ularly  shaped  islands.  In  an  IOPS  application,  REFDIF 
is  employed  to  produce  linear  transfer  functions 
for  refraction  and  shoaling  that  are  used  to  modify  off¬ 
shore  (WAM)  directional  wave  spectrum  to  wave  con¬ 
ditions  found  very  nearshore,  suitable  for  initializing 
the  surf  model. 

Running  wave  models  such  as  REFDIF  each  time 
for  new  boundary  forcing  can  be  a  time  consuming 
process.  A  transfer  function  approach  can  be  applied 
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for  wave  hindcasts,  nowcasts,  or  forecasts  and  can 
greatly  reduce  run  time.  In  this  approach,  calculations 
are  made  in  advance  for  all  possible  frequency  and 
angular  components,  and  the  results  are  saved  in  a  tab¬ 
ular  form.  A  set  of  transfer  functions  (one  each  for 
shoaling  and  refraction)  are  produced  for  a  given  loca¬ 
tion,  and  remain  valid  until  significant  changes  in  the 
nearshore  depth  occur,  for  example,  after  a  major  storm. 
The  transfer  function  consists  of  two  rectangular  matri¬ 
ces.  Each  element  of  the  matrices  corresponds  to  a 
unique  wave  frequency  and  direction.  The  first  matrix 
consists  of  amplitude  coefficients,  representing  the 
change  in  amplitude  for  a  wave  with  a  given  frequency 
and  direction  from  the  offshore  boundary  of  the  model 
to  a  nearshore  location.  The  second  matrix  consists  of 
directions  representing  the  change  in  direction  for  a 
given  wave.  To  modify  a  complete  directional  wave 
spectrum,  the  amplitudes  and  directions  of  individual 
wave  components  are  modified  by  the  corresponding 
transfer  functions  and  then  linearly  combined  with 
results  from  other  components.  An  example  of  this 
technique  is  shown  in  Figure  5  in  which  an  incoming 
wave  at  a  water  depth  of  100  m  with  a  period  of  10  sec¬ 
onds  and  an  offshore  wave  direction  of  -30°,  is  refract¬ 
ed  to  a  wave  angle  of  -14.3°  at  a  water  depth  of  9  m. 


O'Reilly  and  Guza  (1993)  have  successfully  applied 
this  superposition  approximation  (valid  only  where 
there  is  no  wave  breaking)  to  provide  swell  forecasting 
in  the  Southern  California  Bight.  Kaihatu  et  al.  (1998) 
examined  phase-resolving  models  for  a  narrow  shelf 
and  smooth  bathymetry  near  Oceanside,  California, 
and  an  East  Coast  site  at  Camp  Lejeune,  North 
Carolina,  with  fairly  complex  bathymetry  and  distinct 
bottom  features.  They  concluded  that  phase-resolving 
models  such  as  REFDIF  offer  an  accurate  description  of 
monochromatic  swell  wave  propagation  over  highly 
variable  bathymetry,  but  care  must  be  given  in  how  the 
incoming  wave  spectra  are  discretized.  Since  REFDIF 
does  not  account  for  wind-induced  processes,  the 
model  is  considered  valid  only  during  periods  when 
wind  flows  are  of  moderate  strength  or  less,  or  in  a 
fetch-limited  situation,  where  the  model  is  applied  to  a 
relatively  small  geographic  area  (e.g.  a  coast  with  a  nar¬ 
row  continental  shelf). 

Tide  Models 

A  typical  IOPS  implementation  utilizes  a  tide 
model  to  provide  tidal  predictions  for  surf  locations  of 
interest.  IOPS  provides  two  tide  model  options: 
ADCIRC  or  PCTIDES.  ADCIRC  is  run  on  an  unstruc- 
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Figure  5.  Example  of  REFDIF  refraction  coefficients;  incoming  wave  has  a  period  of  10  seconds  and  wave  direction  of  -30°; 
the  modified  wave  direction  is  -14.3°. 
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tured  grid,  requiring  more  expertise  by  the  user  to 
properly  set  up  the  grid  and  apply  appropriate  tidal 
constituent  boundary  conditions.  An  advantage  in  run¬ 
ning  ADCIRC,  however,  is  that  its  triangular  mesh 
allows  one  to  resolve  coastal  features  and  inlets  at  very 
high  resolutions  (on  the  order  of  meters,  if  desired), 
while  a  region  far  removed  from  the  area  of  interest  can 
be  resolved  on  scales  of  tens  of  kilometers.  PCTIDES  is 
much  easier  to  set  up  since  the  grid  spacing  is  equidis¬ 
tant  (Ax=Ay).  The  following  section  provides  back¬ 
ground  information  on  these  models. 

ADCIRC 

ADCIRC-2DDI  (Luettich  et  al.,  1992)  is  a  2-D 
barotropic  (depth-averaged)  hydrodynamic  model  that 
calculates  water  surface  elevations  and  depth-aver¬ 
aged  currents  in  the  nearshore  on  a  finite-element  grid. 
The  inputs  required  to  run  ADCIRC  include: 

♦  Bathymetry  and  shoreline  position 

♦  Specification  of  a  finite-element  grid  and 
boundary  conditions, 

♦  Tidal  constituents  along  the  grid  boundary, 
and 

♦  Wind  speed  and  direction  and  atmospheric 
pressure. 

ADCIRC  has  the  capability  of  activating  the  wet¬ 
ting  and  drying  of  computational  grid  cells.  The  model 
is  also  used  for  both  tidal  and  storm  surge  simulations 
(Blain  et  al.,  1994).  A  more  detailed  description  of 
ADCIRC  can  be  found  in  this  publication  (Blain  et  al., 
this  issue). 

In  an  IOPS  framework,  when  a  surf  model  run  is 
being  set  up,  the  ROAMER  relational  database  is 
queried  to  determine  if  an  ADCIRC  grid  supports  the 
set  of  water  level  locations  requested  by  the  user.  If  an 
ADCIRC  grid  exists,  a  series  of  locations  (latitude,  lon¬ 
gitude)  are  passed  to  an  ADCIRC  post-processor, 
which  interpolates  ADCIRC  water  elevations  from  the 
finite-element  grid  to  the  desired  location.  If  the 
requested  grid  points  do  not  fall  within  the  ADCIRC 
domain,  then  another  tidal-prediction  model,  such  as 
PCTIDES  is  used. 

PCTIDES 

PCTIDES  is  a  globally  relocatable  tidal-prediction 
model,  which  consists  of  a  2-D  barotropic  ocean  model. 
Originally  designed  for  a  PC-based  environment, 
PCTIDES  has  been  configured  to  run  on  UNIX  plat¬ 
forms.  Surface  winds  and  pressures  and/or  astronom¬ 
ical  tides  are  used  to  force  the  model.  A  global  tide 
model,  the  Finite  Element  Solutions  95.1/2.1 
(FES95. 1/2.1),  is  used  to  provide  tidal  conditions  on 
open  boundaries.  All  databases,  except  for  wind  and 
pressure  forcing,  are  internal  to  PCTIDES.  These  data¬ 
bases  include:  (1)  bathymetry,  a  3-minute  interpolated 
version  of  DBDB-V,  (2)  the  FES95. 1/2.1  solutions,  and 
(3)  a  tidal  station  database  consisting  of  more  than  4200 
stations  from  the  International  Hydrographic  Office 
(IHO)  database.  The  IHO  data  can  be  used  for  either 


model  validation  or  for  data  assimilation. 

PCTIDES  can  be  configured  to  run  " stand-alone" 
or  in  a  nested  fashion.  Generally,  a  nested  region  covers 
an  area  smaller  than  the  host  grid  but  with  a  resolution 
3  to  4  times  higher.  The  IOPS  user  can  select  wind  forc¬ 
ing  from  NOGAPS,  which  provides  truly  global  cover¬ 
age  with  a  horizontal  resolution  of  1.0°.  NOGAPS  is 
used  if  wind  forcing  is  important,  but  no  higher-resolu¬ 
tion  wind  fields  are  available.  However,  the  U.S.  Navy 
provides  higher-resolution  COAMPS  forcing  fields  for 
many  geographic  regions  of  interest.  More  detailed 
information  about  PCTIDES  can  be  found  in  this  issue 
(see  Blain  et.  al). 

Under  IOPS,  PCTIDES  is  used  to  provide  a  water- 
level  time  series  (e.g.  a  48-hour  forecast)  for  any  loca¬ 
tion  where  a  surf  forecast  is  to  be  made. 

Surf  Modeling 

The  NSSM  has  been  used  extensively  throughout  the 
U.S.  Fleet  since  it  was  first  introduced  in  the  late  1980s. 
The  model  is  the  primary  software  for  Navy  operational 
surf  forecasting,  and  it  is  also  used  for  system  develop¬ 
ment  and  surf  climate  descriptions.  NSSM  has  gone 
through  several  improvements.  Earle  (1999)  provides  a 
model  overview  and  examples  of  its  applications. 
Further  application  of  NSSM  coupled  to  other  wave  and 
tide  models  can  be  found  in  Allard  et  al.  (1999). 

NSSM  produces  outputs  that  meet  the  operational 
requirements  as  described  in  the  Joint  Surf  Manual 
(Commander,  Naval  Surface  Force,  Pacific  and 
Commander,  Naval  Surface  Force,  Atlantic,1987). 
Model  outputs  include  the  following: 

♦  A  forecast  summary,  which  lists  the  surf  zone 
width,  the  maximum  wave  height  and  long 
shore  current,  direction  of  breakers,  breaker 
type,  breaker  period,  and  the  Modified  Surf 
Index  (MSI),  which  is  a  dimensionless  number 
that  characterizes  overall  surf  conditions. 

♦  An  optional  listing  of  depth,  wave  height, 
wavelength,  percent  of  breaking  waves,  and 
longshore  current  as  function  of  distance  off 
shore. 

Computations  of  Wave  Height  and 
Longshore  Current 

Hsu  et  al.  (2000)  present  an  updated  description  of 
the  NSSM  theory,  equations,  and  numerical  methods; 
thus,  only  a  brief  description  of  the  model  is  given  here¬ 
in.  The  model  is  parametric  and  one-dimensional,  and 
is  largely  based  on  Thornton  and  Guza  formulation 
(1983,  1986).  The  wave  height  in  the  surf  zone  is  com¬ 
puted  based  on  the  energy  flux  balance.  The  energy 
flux  equation  balances  the  wave  and  roller  fluxes  of 
energy  density  to  energy  dissipation  due  to  the  depth- 
induced  breaking.  The  model  uses  wave-roller  dissipa¬ 
tion  developed  by  Lippmann  et  al.  (1996).  Main 
assumptions  in  the  surf  model  are:  (1)  straight  and  par¬ 
allel  bottom  contours  and  (2)  narrow  banded  waves 
(i.e.  single-peak  wave  frequency). 
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Longshore  current  calculations  are  based  on  radia¬ 
tion-stress  longshore  current  theory  first  developed  by 
Longuet-Higgins  (1970a,b).  The  momentum  conserva¬ 
tion  equation  considers  the  production  of  momentum 
fluxes  from  breaking  waves  and  wind,  horizontal  tur¬ 
bulent  mixing,  and  bottom  friction.  Recently,  Hsu  et  al. 
(2000)  have  shown  that  the  longshore  current  estimates 
can  be  considerably  improved  by  making  relatively 
simple  changes  to  bottom  friction  formulation.  After 
the  bottom  friction  coefficient  in  the  model  is  changed 
to  a  depth-dependent  bottom  friction  function,  the 
longshore  current  estimates  are  significantly 
improved.  It  should  be  noted  that  the  variable  bottom 
friction  function  reflects  the  shoreward  increase  in 
friction  due  to  sediment  sorting  and  compensates 
for  the  lack  of  vertical  eddy  diffusivity  in  one-dimen¬ 
sional  models. 

NSSM  Model  Validation 

NSSM  has  been  validated  extensively  (Mettlach  et 
al,  2000)  with  data  collected  from  field  studies  and  lab¬ 
oratory  experiments.  Figure  6a  shows  a  plot  of  long¬ 
shore  current  from  measurements  taken  during  the 
DELILAH  and  DUCK94  experiments  versus  NSSM 
using  optimal  bottom  friction  and  horizontal  mixing 
coefficients.  Figure  6b  shows  very  good  agreement  with 
NSSM  wave  height  versus  observations.  This  study 
shows  that  with  accurate  inputs  to  the  model,  NSSM 
can  provide  skill  in  the  depiction  of  surf  conditions. 


NSSM  Inputs 

IOPS  is  an  integrated  system  providing  the 
hooks  required  to  generate  a  surf  forecast.  The  follow¬ 
ing  section  describes  the  necessary  inputs  required 
to  run  NSSM. 

Directional  Wave  Spectra:  NSSM  is  initialized 
with  directional  wave  spectra  obtained  from  one  of  the 
two  shallow-water  wave  models,  STWAVE  or  REFDIF. 
Spectra  are  saved  at  water  depths  of  approximately  7-9 
m,  just  outside  of  the  surf-zone  where  wave  breaking 
has  not  occurred. 

Tides:  Water  level  information  provided  by 
ADCIRC  or  PCTIDES  is  input  to  NSSM  to  adjust  the 
beach  profile.  Surf  forecasts  are  generally  made  for  48- 
hour  periods;  a  tidal  prediction  with  a  temporal  resolu¬ 
tion  of  at  least  6  hours  is  most  commonly  used. 

Winds:  10-m  surface  winds  fields  are  input  into 
NSSM  and  used  in  the  calculation  for  the  longshore 
current.  The  default  model  for  wind  input  is  NOGAPS. 
Otherwise,  CO  AMPS  is  used. 

Beach  Profiles:  The  user  has  several  options  avail¬ 
able  to  incorporate  a  beach  profile(s)  into  the  model.  A 
nearshore  depth  profile  extends  from  land  to  a  water 
depth  of  about  8  m.  Some  of  the  beach  profile  options 
are  discussed  below: 


Measured  Longshore  Current  (meters/sec)  Measured  Wave  Height  (meters) 

Figure  6.  Scatter  diagram  of  a)  measured  longshore  current  vs.  corresponding  model  estimates  of  current.  Observed  data  are 
from  DELILAH  and  DUCK94  measurements.  Based  on  1930  observations ,  the  linear  correlation  coefficient  is  0.89  with  a 
standard  deviation  of  0.22  m/s.  b)  measured  root-mean  wave  heights  vs.  estimates  of  wave  height.  Based  on  1796  observations , 
the  linear  correlation  coefficient  is  0.945  with  a  standard  deviation  of  0.15  m. 
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Equilibrium  Sediment  Profiles 


Figure  7.  Equilibrium  beach  profiles 


Equilibrium  Beach  Profile:  An  equilibrium  beach 
profile  (Dean,  1977)  based  on  sediment  type  (and  grain 
size)  is  the  default,  if  no  other  information  is  available. 
In  the  Dean  study,  502  beach  profiles  from  the  U.S.  East 
Coast  and  Gulf  of  Mexico  were  examined.  Figure  7 
shows  equilibrium  beach  profiles  based  on  grain  sizes 
ranging  from  pebbles  to  very  fine  sand  and  clay. 


Imported  Beach  Profiles:  IOPS 
allows  a  user  to  import  beach  profiles 
collected  from  a  variety  of  sources 
including  Sea  Air  Land  (SEAL)  team 
beach  surveys  and  those  collected 
from  systems  such  as  the  Scanning 
Hydrographic  Operational  Airborne 
Lidar  Survey  (SHOALS).  During  .the 
Rapid  Environmental  Assessment 
phase  of  NATO  Linked  Seas  2000 
Exercise  near  Pinheiro  da  Cruz, 
Portugal,  SHOALS  was  used  to  col¬ 
lect  beach  profiles.  Operating  on  a 
fixed-wing  Twin  Otter  aircraft  at  an 
altitude  of  300  m,  SHOALS  (Lillycroft 
et  al.,  2000)  was  used  to  collect 
nearshore  bathymetry  with  a  resolu¬ 
tion  for  each  beach  sounding  of  4  m. 
The  NAVOCEANO  WSC  utilized  the 
SHOALS  data  in  providing  wave  and 
surf  forecasts  to  support  the  exercise. 

NSSM  Output 

An  example  of  NSSM  output 
depicted  in  Figure  8  a-e  shows  a  48- 
hour  time  series  of  surf-related 
parameters  from  the  IOPS  coupled 
wave-tide-surf  modeling  system.  If  a 
landing  is  being  considered  for  several  beaches,  this 
forecast  can  be  used  as  a  tool  to  determine  the  best  loca¬ 
tion  based  on  predicted  surf  conditions.  A  text  version 
of  NSSM  output  is  shown  in  Table  1.  The  Modified  Surf 
Index  (see  Table  1)  can  be  used  by  mission  planners  to 
determine  if  surf  conditions  are  near  or  exceed  critical 
levels  based  on  the  type  of  landing  craft. 

Summary  and  Future  Plans 

An  integrated  modeling  approach  has 
been  described  in  which  the  deep-water 
WAM  is  coupled  to  shallow-water  wave, 
tide,  and  surf  models.  The  system  is 
designed  to  be  modular  in  fashion  so  that 
improvements  or  replacements  for  the 
modules  can  be  incorporated.  IOPS  soft¬ 
ware  and  documentation  was  delivered  to 
NAVOCEANO  in  December  2001.  Testing 
and  evaluation  of  the  system  began  in 
January  2002. 

Future  plans  include  integrating  the 
Simulating  WAves  Nearshore  (SWAN) 
wave  model  into  IOPS.  SWAN  includes 
improved  approaches  for  wave  propaga¬ 
tion  and  an  ability  to  operate  efficiently  on 
very  fine  grid  meshes  (order  100  m  or  less). 
It  is  a  full  plane  model  (onshore  and  off¬ 
shore  winds  and  waves)  and  has  both 
time-dependent  and  steady-state  modes. 
The  model  is  designed  to  nest  with  WAM. 
SWAN  has  undergone  significant  testing 


Table  1 

Standard  NSSM  output 


Coded  Surf  Forecast  Follows 


Significant  Breaker  Height 

alfa  =  1 .46  m 

Maximum  Breaker  Height 

bravo  =  2.26  m 

Dominant  Breaker  Period 

charlie  =  1 0.5  sec. 

Dominant  Breaker  Type 

delta  =  Spilling  Surf 

(56%  Spilling,  44%  Plunging,  0%  Surging) 

Breaker  Angle  (toward  left  flank) 

echo  =  1 .9  deg. 

Littoral  Current  (toward  left  flank) 

foxtrot  =  0.2  m/s 

Number  of  Surf  Lines 

golfl  =  1 .8 

Surf  Zone  Width 

golf2  =  74.08  m 

Wind  Speed 

hotel  1  =  5  m/s 

Wind  Direction 

hotel2  =  1 26.4  deg. 

Modified  surf  index  =  3.5 
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Figure  8.  Surf  forecast  for  a  potential  landing  beach:  a)  beach  profiles  used  in  the  surf  calculation,  b)  tidal  prediction  with  fore¬ 
cast  10-m  winds,  c)  significant  /maximum  breaker  heights  within  the  surf  zone,  d)  surf  zone  width  and  breaker  type,  e)  long¬ 
shore  current  and  f)  breaker  angle. 


(perhaps  the  most  extensive  testing  of  any  shallow- 
water  wave  model)  and  has  an  international  user  base 
of  over  100  users.  The  model  is  well  documented. 
SWAN  is  being  run  to  support  the  Northern  Gulf 
Littoral  Initiative  (NGLI)  in  the  Mississippi  Bight. 
Plans  are  underway  to  begin  an  evaluation  and  com¬ 
mence  testing  of  SWAN  at  NAVOCEANO. 

A  distributed  IOPS  (DIOPS)  has  been  installed  at 
the  Navy  Pacific  Meteorology  Oceanography  Center  in 
San  Diego,  California.  DIOPS  allows  a  user  to  run  a 
host  of  models  in  an  object-oriented  environment  on  a 
UNIX  workstation.  Atmospheric  forcing  and  bathyme¬ 
try  are  accessed  through  the  U.S.  Navy's  Tactical 
Environmental  Data  Server  (TEDS).  It  is  anticipated 
that  directional  wave  spectra  from  NAVOCEANO's 
Regional  WAM  model  will  be  distributed  via  TEDS  to 
initialize  DIOPS  shallow-water  wave  models,  ultimate¬ 
ly  providing  situational  awareness  of  the  littoral  and 
surf  environmental  conditions  to  the  warfighter  on- 
demand.  S3 
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